Hall effect of YbRh 2 Si2 and relatives in the light of electronic structure calculations 
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We report experimental and theoretical investigations of the Hall effect in YbRli2Si2 and its ref- 
erence compounds LuRli2Si2 and YbIr2Si2- Based on band-structure calculations we identify two 
bands dominating the Hall coefficient in all these compounds. For the case of LuRli2Si2 — the non- 
magnetic reference compound of YbRli2Si2 — the temperature dependence of the Hall coefficient is 
described quantitatively to arise from two hole-like bands. For YbIr2Si2 and YbRli2Si2, renormalized 
band calculations yield two bands of opposite character. In YbRli2Si2 these two bands almost com- 
pensate each other. We present strong indications that sample dependences of the low-temperature 
Hall coefficient observed for YbRli2Si2 arise from slight variations of the relative scattering rates of 
the two bands. Minute changes of the composition appear to be the origin. The results of our band 
structure calculations reveal that a transition of the 4/ electrons from localized to itinerant leads 
to a decrease of the Hall coefficient. 

PACS numbers: 71.27. +a, 71.18.+y, 72.20.My, 73.43.Nq 
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I. INTRODUCTION 

The heavy-fermion compound YbRh2Si2 has emerged 
as a prototypical system for the investigation of quan- 
tum critical phenomena.- Pronounced non- Fermi-liquid 
properties arise due to the proximity to a quantum crit- 
ical point (QCP). 2 In its ground state, YbRli2Si2 or- 
ders antiferromagnetically below the Neel temperature, 
Tn = 70mK. 3 By applying a small magnetic field of 
B c = 60 mT within the basal plane, the magnetic or- 
der is suppressed to zero temperature, thus accessing the 
field-induced QCP4. 

Hall-effect measurements turned out to be of central 
importance to understand the nature of the QCP as they 
allow to discriminate two different theoretical scenarios 5 
as discussed below. The Hall coefficient i?n of YbRh2Si2 
was measured as the compound was driven from the mag- 
netically ordered state across the QCP towards the LFL 
regime by increasing the magnetic field Since anoma- 
lous contributions^ are negligible^ at low temperatures, 
i?H is directly related to the Fermi surface volume. The 
Hall coefficient was found to exhibit a crossover linked to 
the QCP which resides on top of a smooth background. 
Since this crossover sharpens to a discontinuous jump in 
the extrapolation to zero temperature these results imply 
an abrupt change of the Fermi surface at the QCP. Such a 
Fermi surface reconstruction is at variance with the pre- 
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dictions of the standard spin-density-wave theoryji£~— 
Rather, the results suggest a new class of theoretical de- 
scriptions to be applied in YbRli2Si2, namely the Kondo- 
breakdown scenario in which the 4/ electrons are itiner- 
ant on the high-field side of the QCP only^^-i 5 - Con- 
sequently, the Hall effect represents a key experiment to 
identify the unconventional nature of the quantum criti- 
cality in YbRli2Si2. 

On the other hand, it was pointed out that the Hall 
coefficient of YbRli2Si2 is not simply proportional to the 
inverse charge carrier concentration since the assumption 
of a spherical Fermi surface with a single band at the 
Fermi energy E-p is not valid in this material as shown 
by various band-structure calculation a 16 i 17 and photoe- 
mission studiesi 16 ' 18 Several calculations yield multiple 
bands crossing E-p with canceling positive and negative 
contributions to the Hall coefficientj 17 i 19 Thus, the cal- 
culated Hall coefficient critically depends on the method 
used, and it remains an outstanding challenge to inter- 
pret the measured Hall coefficient quantitatively in terms 
of band-structure calculations. 

The Hall-effect measurements 2 ^ performed in a low- 
temperature setup with improved resolution on a num- 
ber of YbRli2Si2 single crystals reproduced the results 
on a crystal used in Ref. la. Other crystals of different 
quality, however, show strong sample dependences be- 
low 20 K. Although the critical crossover is found to be 
virtually independent of sample quality, it remains to be 
understood why the underlying background exhibits such 
strong sample dependences^ 

Here, we present experimental and theoretical progress 
which helps to refine our understanding of the Hall effect 
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in YbRti2Si2. The issue of the low-temperature sam- 
ple dependences as well as the characteristics of R^(T) 
in YbRli2Si2 are addressed. Comparison with the Hall- 
effect data of both the non-magnetic reference compound 
LuRli2Si2 and the heavy- fermion compound YbIr2Si2 al- 
lows us to discriminate various contributions to the Hall 
coefficient. Our renormalized band-structure calculations 
yield excellent agreement with the experimentally deter- 
mined Hall coefficient for YbIr2Si2 and LuRli2Si2. They 
provide a reliable basis to understand the sample depen- 
dences in YbRh 2 Si2- Moreover, enable us to relate the 
field induced crossover of the Hall coefficient to a change 
in the carrier concentration^. 

We present the details and the results of the electronic 
structure calculations in section[H] This includes the cal- 
culation of the Hall coefficient in section III El In section 
IIIII the results of the electronic transport measurements 
on LuRli2Si2, YbIr2Si2, and YbRli2Si2 are presented and 
discussed in the light of determined electronic band struc- 
tures. 



II. BAND-STRUCTURE CALCULATIONS 
A. Models of electronic structure 

1. Local moment regime 

We begin by calculating the electronic structure of 
YbRli2Si2 and YbIr2Si2, assuming that the Yb ions are 
in the 4/ 13 configuration. With this approximation we 
model the Fermi surface and the quasiparticle bands in 
the local moment regime. As there are exactly 13 4/ 
electrons or one 4/ hole per Yb site the single-particle 
excitations of the 4/ shell involve valence transitions 
4/ 13 -)• 4/ 14 and 4/ 13 -> 4/ 12 which occur at high 
energies only. Consequently, the 4/ degrees of freedom 
do not contribute to the low-energy excitations in the 
vicinity of the Fermi surface. In this energy range, the 
single-particle excitations are derived from the weakly 
correlated (non-/) conduction states which form coher- 
ent Bloch states. We determine the dispersion of these 
bands by standard band-structure calculations. The ef- 
fective potentials are generated self-consistently within 
the Local Density Approximation (LDA) to density func- 
tional theory. The strong Coulomb repulsion among the 
4/ electrons which suppresses charge fluctuations is im- 
plicitly accounted for by treating the 4/ electrons as part 
of the ion core assuming that they do not hybridize with 
the conduction states. This assumption seems justified 
for the systems under consideration whose 4/ valence de- 
viates only weakly from the integer value. We refer to this 
method as /-core calculation. By using the /-core calcu- 
lation for YbRli2Si2 to interpret the results obtained on 
LuRli2Si2 we rely on the facts that the lattice parame- 
ters agree within the experimental error— and that the 
results of the /-core calculation are independent of the 
4/-occupancy. 



The partially filled / shell of the 4/ 13 -configuration 
necessarily carries a magnetic moment in agreement with 
Kramers' theorem. The presence of local magnetic mo- 
ments is reflected in the Curie- Weiss behavior observed 
at elevated temperatures in the magnetic susceptibility 
of YbRli2Si2 and YbIr2Si2- The 4/ moments, however, 
interact only weakly with the conduction states as can 
be inferred, e.g.. from the low magnetic ordering temper- 
ature in YbRh.2Si2. We neglect the potential reconstruc- 
tion of the conduction electron Fermi surface that may 
result from the long-range antiferromagnetic order and 
account only for the 4/ charge which contributes to the 
potential seen by the conduction electrons. This amounts 
to effectively averaging over the local magnetic degrees 
of freedom in determining the self-consistent potentials. 
When comparing with experiment the bare bands de- 
rived from the effective potentials have to be renormal- 
ized by local 4/ excitations. Scattering off Crystalline 
Electric Field (CEF) excitations may enhance the effec- 
tive masses and reduce the life-times of the conduction 
electrons. With these effects properly accounted fo r 22 i 23 
the /-core model should quantitatively describe the elec- 
tronic properties of Yb-based heavy-fermion compounds 
at elevated temperatures. However, for the Hall effect 
of the heavy-fermion compounds one has to take anoma- 
lous contributions into account which arise from the skew 
scattering at the local /-moments. These contributions 
may only be neglected at very low temperatures where, 
on the other hand, the /-core calculation is insufficient 
for the description of heavy-fermion compounds. Rather, 
we shall use the Renormalized Band Calculation (RBC) 
to understand the Hall coefficient in the heavy-fermion 
compounds. 

Treating the 4/ electrons as part of the ion core can 
be viewed as an extreme limit of an LDA+U calculation. 
Therefore, we shall compare our data with recent results 
obtained from LDA+U (Ref. [fi). The LDA+U calcula- 
tion explicitly includes the magnetic moments of the 4/ 13 
configuration assuming long-range ferromagnetic order. 
This treatment preserves the translational invariance of 
the underlying lattice. It removes, however, the spin de- 
generacy of the conduction bands as they are split by the 
Zeeman effect. This splitting is rather small reflecting 
the weak coupling between the 4/ states and the conduc- 
tion electrons. For this reason, we anticipate the energy 
bands of the LDA+U and the /-core calculation to agree 
in the low-energy regime, i.e., in the vicinity of the Fermi 
surface. 



2. Heavy Fermi liquid regime 

The strongly renormalized heavy quasiparticle bands 
are determined by means of the renormalized band 
metho d 2 ' 25 which combines material-specific ab initio 
methods and phenomenological considerations in the 
spirit of Landau. The key idea is to construct an ef- 
fective Hamiltonian for the low-energy excitations which 
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uses the ab initio potentials for the weakly correlated con- 
duction electron channels while introducing one param- 
eter to account for the specific local correlations among 
the 4/ electrons. The parameter is determined once by 
fitting to experiment and is kept fixed during subsequent 
investigations. A detailed description of the method and 
typical results for Ce-based compounds are given in Refs. 
[26l . I27I . Operationally, it amounts to transforming the 
/-states of the spin-orbit ground state multiplet at the 
lanthanide site into the basis of CEF eigenstates \m) and 
introducing resonance-type phase shifts 



fjf(E) ~ arctan \ (1) 
E-e f 

where the resonance width A f accounts for the renormal- 
ized quasiparticle mass. The resonance energies et m = 
£/ + &m refer to the centers of gravity of the /-derived 
quasiparticle bands. Here if denotes the position of the 
band center corresponding to the CEF ground state while 
5 m are the measured CEF excitation energies. One of 
the remaining two parameters, e/, is determined by im- 
posing the condition that the charge distribution is not 
altered significantly by introducing the renormalization. 
This makes the RBC a single-parameter scheme. The 
free parameter, Af. is adjusted in such a manner that 
the coefficient of the linear-in-T specific heat at low tem- 
peratures is reproduced. The effective band structure 
Hamiltonian constructed along these lines corresponds 
to a hybridization model which closely parallels the one 
obtained from the periodic Anderson model in mean-field 
approximation. Alternatively, the RBC can be viewed as 
a parametrization scheme for the variation with energy 
of the real part of the local 4/ electron self-energy. The 
parameter to be determined by experiment is the slope 
at the Fermi energy while the value at E-p is fixed by 
retaining the charge distribution. 

The method has been shown to reproduce Fermi sur- 
faces and anisotropics in the effective masses of a great 
variety of Ce-based compounds. In addition, it allows to 
predict Fermi liquid instabilities^— . 

In calculating the coherent 4/-derived quasiparticle 
bands in Yb-based heavy-fermion compounds we essen- 
tially follow the procedure for the Ce case as described 
above. We have to account for the fact that Yb can be 
considered as the hole analogue of Ce. Operationally this 
implies that we have to renormalize the 4/ j = 7/2 chan- 
nels at the Yb sites instead of the 4/ j = 5/2 states in the 
Ce case. As the 4/ hole count is slightly less than unity 
the center of gravity e / will lie below the Fermi energy. 
In addition, we have to reverse the hierarchy of the CEF 
scheme, i.e., 

e f < ; i fm = e f - S m . (2) 



B. Computational method 

The calculations are done on the basis of the experi- 
mental lattice parameters a = b = 4.007 A, c = 9.858 A 
for YbRli2Si2 and LuRli2Si2(cf. section III A 1[) and a — 
b = 4.032 A, c = 9.826 A for YbIr 2 Si 2 (I-type)^ The 
band structures were obtained by the fully relativistic 
formulation of the linear muffin-tin orbitals method^ir— . 
We adopt the atomic-sphere approximation including the 
combined correction term which contains the leading 
corrections 3 ^. In solving the band-structure problem, we 
include s-p-d-/-components at the Yb and the transi- 
tion metal (Rh, Ir) sites and s-p-c?-components at the 
Si sites. The spin-orbit interaction is fully taken into 
account by solving the Dirac equation. Although the 
relativistic effects hardly change the electron density dis- 
tribution they nevertheless influence the actual location 
of the energy bands. This aspect is particularly impor- 
tant for the renormalized band structure since the spin- 
orbit splitting of the c?-states is rather large on the en- 
ergy scales relevant for the strongly renormalized heavy 
quasiparticles. Exchange and correlation effects were in- 
troduced using the Barth-Hedin potential^!. The band 
structure was converged for 405 k-points within the irre- 
ducible wedge, whose volume equals 1/16 of the Brillouin 
zone. The density of states (DOS) was evaluated by the 
tetrahedron method with linear interpolation for the en- 
ergies. For the conduction band the DOS was calculated 
at 0.25 mRy (« 3.4 meV) intervals. To obtain reliable 
values for the transport integrals the energies were cal- 
culated at 2601 k-points within the irreducible wedge. 
Subsequently, the bands were interpolated using Mathc- 
matica VI, and the result was used to numerically eval- 
uate the desired quantities. 



C. Electronic structure 

1. YbRh,2Si2 in the local-moment regime 

Figure [1] displays the electron bands of YbRh.2Si2 in 
the vicinity of the Fermi energy along symmetry lines 
with the 4/ electrons being treated as part of the ion 
core. Here, the band states have predominantly Rh 4c? 
character with some admixture of Yb 5c? character. The 
dispersion of YbRli2Si2 agrees rather well with the re- 
sults of recent LDA+U calculation s 16 ' 39 . In addition, it 
is consistent with energy bands deduced from photoemis- 
sion studies J£ 

In the 4/-core calculation, the broad bands intersecting 
the Fermi energy are exclusively formed by the non-/ 
conduction states. This is reflected in the low DOS at 
the Fermi energy N(E-p) — 2.1 states/ (eV unit cell) for 
the /-core calculation of YbRli2Si2 as shown in Fig. [5J 
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FIG. 1: YbRh2Si2: Electronic bands along symmetry lines 
with the Fermi energy Ef = chosen as reference energy. The 
Yb 4/ electrons are treated as part of the ion core. We follow 
the notation of Ref. [H using the labels Z (0,0,1), F (0,0,0), 
X (1,1,0), P (1,1,1) and N (1, §, §) in units of (=,_=,=). 
The labels s, f, and u refer to (5,0,0) , (6,0,2) and (6,6,2), 
respectively, with 5 = 1 + (- ) 2 and 6 = 1- (- ) 2 . 
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FIG. 2: (color online) YbRh2Si2: Comparison of the total 
DOS in the local moment regime (/-core calculation, solid 
line) and in the heavy- Fermi-liquid regime (RBC, red shaded 
area). The reference energy is the Fermi energy Ef = 0. The 
two bands in the low-energy part are derived from the Si s 
states. The dominant features are the Rh Ad bands which 
hybridize with Si p and Rh s states near the bottom of the d 
bands and with Si d and Rhp states near the top, respectively. 



state which is a superposition of \ j = 7/2; j z = ±5/2) and 
\j = 7/2; j z = +3/2) and which is well separated from the 
excited states. The CEF parameters and the CEF eigen- 
states are given in Ref. |42|. Using the effective quasipar- 
ticle resonance widths of A/ ~ 20K and Af ~ 40K as in- 
ferred from specific heat and thermopower measurements 
for the Rh (Ref. HHH) and Ir (Ref. [Hi) compounds yields 
the band structures displayed in Fig. [3] The dispersion of 
the renormalized bands of the Rh- and the Ir-compound 
are rather similar, the band widths scale with the charac- 
teristic temperatures. We shall concentrate on the results 
for the Rh-compound in the subsequent discussion. 

The RBC yields narrow /-derived quasiparticle bands 
in the vicinity of the Fermi energy, whereas the disper- 
sion of the non-/ bands is essentially unaffected. This 
can be seen from Fig. [2] in which the DOS derived from 
the renormalized bands are compared with the /-core 
counterpart. The expanded view of the RBC DOS in 
the low-energy regime as depicted in Fig. @] shows the 
contributions of the CEF-split 4/ states. The CEF ex- 
citations appear in the occupied part of the spectrum 
below the Fermi energy. The hybridization and hence 
the effective quasiparticle masses are rather anisotropic. 
The renormalized band calculations yield a DOS of 
290 states/ (eV unit cell) at Ep corresponding to specific 
heat coefficient eSOmJmoF 1 K~ 2 . For YbIr 2 Si 2 a DOS 
of 48 states/ (eV unit cell) at Ep is calculated correspond- 
ing to a Sommerfeld coefficient of HSmJmoP 1 K~ 2 . 



D. Fermi surface and quasiparticles 

1. Local-moment regime 

The LDA calculation for localized 4/ electrons pre- 
dicts three bands to cross the Fermi energy and leads 
to the Fermi surface which closely resembles previous 
results 15 - 44 . It consists of three separate sheets. The 
two main sheets form a hole surface centered around the 
Z point, and a complex, multi-connected surface. Fol- 
lowing Ref. [16| we shall refer to them as 'donut' and 'jun- 
gle gym', respectively. In addition, there is a small T- 
centered electron surface, the 'pill box'. We shall focus 
on the two main sheets which are displayed in Fig. [5] as 
these two dominate the electronic properties. 



2. YbRh2Si2 and YbIr2Si2 in the heavy-Fermi-liquid regime 

The calculations reported here adopt a CEF scheme 
which is consistent with susceptibility and inelastic 
neutron-scattering dat a 40 ' 41 . The latter indicate that 
the 4/ 13 states in YbIr2Si2 and YbRli2Si2 are split 
into 4 doublets with the energies 0— 18 — 25 — 36 meV 
and — 17 — 25 — 43meV, respectively. The low-energy 
properties are mainly determined by the CEF ground 



2. Heavy-Fermi-liquid regime 

The RBC also predicts two major sheets at the 
Fermi surface whose topologies resemble those found by 
LDAl&i^ (see Fig. [5J. The major sheets of the Rh- 
compound and its Ir-counterpart are rather similar. The 
main difference occurs in the small pockets: The small 
T-centered electron pocket of the Rh-compound is absent 
in the Ir-System where we find a Z-centered hole pocket 
instead. 
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FIG. 3: Dispersion of the renormalized bands along symme- 
try lines for YbRli2Si2 (upper panel) and YbIr2Si2 (lower 
panel). The anisotropy of the CEF ground state leads to 
highly anisotropic hybridization strength which affects the rel- 
ative shifts and the widths of the bands. The topology of the 
Fermi surfaces is mainly determined by the steep conduction 
bands. The symmetry of the CEF ground state is reflected in 
effective mass anisotropies. The coordinates of the symmetry 
points Z, r, X, P, and N are specified in the caption of Fig. [1] 



From the comparison of the /-core results and those 
of the RBC it is obvious that they represent 'small' and 
'large' Fermi surfaces, respectively. The difference of the 
Fermi volume accounts for the additional states related 
to the large quasiparticle DOS at E-p for the heavy-Fermi- 
liquid limit (cf. Fig. |5J|. 



E. Calculation of the Hall coefficient 

For the chosen experimental geometry and using the 
Boltzmann approximation, the Hall coefficient is given 
in the low-field limit byi& 



a xx(i 



(3) 
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FIG. 4: Quasiparticle DOS for YbRli2Si2 in the renormalized 
band structure calculation. For a comparison with the DOS 
plotted in Fig. [1] consider the different scales. 



where the conductivity tensor elements 
<J xx (i) = e 2 ~^~V(i)v£(i,k) ( - 



Of 



dE(i,k) 



(4) 



and 



e 1 

<Jxyz(i) = -nE r2 W [vx(i,k)vy(i,k)M^£(i,k) 



c n 



-^(i,k)M- 1 (*,k)] 



df 



dE(i,k) 



(5) 



are summed over all bands i intersecting the Fermi sur- 
face. We enumerate the 'donut' and the 'jungle gym' with 
i = 1 and i = 2, respectively. We restrict ourself to the 
first two bands which dominate the total conductivities 
and neglect small pockets. In the following, we allow for 
the possibility that the relaxation time r may vary from 
one band to another but we shall neglect the variation of 
the relaxation time r with wave vector k (with compo- 
nents k a ). Here, e and c denote the electron charge and 
the vacuum speed of light, respectively. Q represents the 
volume of the Brillouin zone and / the Fermi distribution 
function. The components of the velocity 



and of the inverse mass tensor 



M^(*,k) = 



1 



<:) 2 



h 2 dk a dki 



-E(i,k) 



(G) 



(7) 



are deduced from the energy bands E(i,k). 

For the discussion of Fermi-surface effects we write the 
longitudinal and transverse conductivity as 



cr X x{i) = cr{i)(Jxx(i) 

<Jxyz{i) = VB{i)Vxyz{i) 



(8) 
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FIG. 5: Calculated Fermi surfaces of YbRli2Si2: Major sheets of the Fermi surface in the /-core calculation representing the 
local- moment regime (top row) and the heavy-fermion regime (bottom row). In addition to the 'donut' (left panel) and the 
'jungle gym' (right panel) there is a small electron surface ('pill box') which is not displayed here. The topology of the Fermi 
surface agrees well with previous findings . 16 ' 19 ' 44 ' 45 



with the prefactors 

e 2 

a(i) = — r(i)n(i) 
m 

a B (i) = K-( T (i)) 2 n(i) (9) 

being the Drude result for a gas of free particles with 
charge |e|. The particle density n(i) corresponds to the 
number of occupied states per unit cell in band i while 
the reduced transport integrals a xx {i) and a xyz {i) ac- 
count for the deviations of the conductivity tensor el- 
ements and the particle density from the free particle 
picture. The results for the systems under consideration 
are summarized in Table Q] 

For the /-core calculations we obtain positive Hall 
(transverse) conductivity for both bands corresponding 
to hole-like character of the charge carriers. The crucial 
point is that for the RBC results, by contrast, the 'jungle- 
gym', is predominantly electron- like as can be inferred 
from the reduced transport integrals listed in Tab. |U 
Moreover, we find for YbRli2Si2 that the two bands al- 
most compensate each other. This is seen by the fact 
that the products n(i)a xyz (i) of the two bands are close 
to each other in magnitude and of opposite sign. Their 



TABLE I: Calculated reduced transport integrals for the two 
different bands (i — 1,2). The results derived for the two dif- 
ferent Fermi-surface models are compared. The Fermi surface 
results for YbIr2Si2 within the RBC are included for compar- 
ison. See text for methods. 

System Method i n(i) a xx (i) a xyz (i) n(i)a xyz (i) 

YbRh 2 Si 2 4/ core 1 1.76 0.197 +0.289 +0.50864 
2 1.22 0.384 +0.153 +0.18666 

YbRh 2 Si 2 4/ RBC 1 1.37 0.0137 +0.00275 +0.0037675 
2 0.63 0.0747 -0.00652 -0.0041076 

YbIr 2 Si 2 4/ RBC 1 1.42 0.051 +0.00323 +0.0045866 
2 0.58 0.138 -0.01003 -0.0058174 



sum determines the numerator of Eq. [3l Since we allow 
for different relaxation rates of the individual bands, this 
gives rise to a weighting of these two terms in the sum of 
Eq. [31 Consequently, the total Hall coefficient very sensi- 
tively depends on the relative relaxation rates of the two 
bands. Even the sign of i?n may change if this balance 
is shifted only slightly toward the electron-like band. We 
shall discuss later that this might relate to the sample 



7 



dependences observed in YbRli2 Si2 . 

The calculated transport integrals vary only slightly 
with the position of the Fermi level. This is contrary to 
the result of Ref. H3 and reflects the different methods 
used: The LDA calculations of Ref. Il7| are not able to 
account for the position of the 4/ level with respect to 
the Fermi energy The RBC on the other hand takes the 
correlation effects into consideration and, thus, does not 
rely on a shift of the 4/ level position. 



III. COMPARISON TO HALL EFFECT 
MEASUREMENTS 

In this section we present the Hall-effect measurements 
and use the above results of the electronic structure cal- 
culations to advance our understanding of the experimen- 
tal observations. 



A. Samples 

Single crystals of YbRIi2Si2, YbIi'2Si2 and LuRli2Si2 
were synthesized applying an In flux-growth technique as 
described earlier' We note that within this work we con- 
centrate on the I- type phase of YbIr2Si2 which is isostruc- 
tural to YbRh 2 Si 2 ^ 

In LuRli2Si2, also isostructural to YbRli2Si2, the Lu 3+ 
has 14 / electrons and consequently retains a fully 
occupied / shell without magnetic moment. There- 
fore, it serves as a non-magnetic reference compound to 
YbRli2Si2. An assignment of the YbRli2Si2 /-core calcu- 
lations to LuRli2Si2 is justified by the fact that LuRli2Si2 
has equal lattice parameters within experimental error. 
This allows us to model the experimentally observed tem- 
perature dependence of the Hall coefficient. The /-core 
calculations yield a DOS (section III C 1|) which corre- 
sponds to a bare linear-in-T specific heat coefficient of 
7 w 5mJmol _ K~ 2 in good agreement with the experi- 
mental value 7 « 6.5mJmol K~ 2 found for LuRli2Si2 
(not shown). The resistivity as displayed in the inset 
of Fig. [HI is approximately linear in T above 100 K with 
p(300 K) = 20 cm. Both the specific heat and the re- 
sistivity indicate that LuRli2Si2 is a simple non- magnetic 
intermetallic compound. 

YbRIi2Si2 and YbIr2Si2 exhibit pronounced heavy- 
fermion behavior in various properties'^ In particular, 
the specific heat is largely enhanced (cf . section III C 21) . 
However, their ground states differ: YbRli2Si2 exhibits 
antifcrromagnetic order at zero magnetic field, whereas 
YbIr2Si2 is paramagnetic obeying Landau Fermi liquid 
(LFL) behavior below 200 mK^ Proximity of YbIr 2 Si 2 
to a QCP is indicated by a logarithmic divergence of 
the specific heat for temperatures above 200 mK. Since 
YbIi - 2Si2 has a larger unit cell volume than YbRIi2Si2 
it is assumed to be located on the paramagnetic side of 
the QCP as unit-cell expansion weakens magnetic order- 
ing in Yb-systems. Consequently, YbIr2Si2 serves as a 



reference compound with fully itinerant 4/-states as ac- 
counted for in the RBC. For YbRli2Si2, in its ground 
state in zero magnetic field, by contrast the / electrons 
appear to be localized as inferred from the Fermi surface 
reconstruction. 6 



B. Experimental setup 

All samples were polished to thin platelets of thick- 
ness in the range 25 pm < i < 80 pm. Subsequently, 
the samples were prescreened via resistivity p(T, B) mea- 
surements to ensure In-free samples. The current / was 
driven within the crystallographic ab plane. The mag- 
netic field B was applied along the c axis, thus, inducing 
the Hall voltage perpendicular to / within the tetrag- 
onal plane, (see inset of Fig. [7]). To measure the Hall 
effect the transverse voltage V y was monitored. In order 
to cancel out magnetoresistance components due to con- 
tact misalignment, the Hall resistivity was obtained from 
the antisymmetric component of the field-reversed trans- 
verse voltage, pn(-B) = t [V y (+B) — V y (-B)] /2I. The 
linear-response Hall coefficient i?n was derived as the 
slope of linear fits to the Hall resistivity pn (B) for fields 
B < 0.4 T. Only the low-temperature Hall resistivity of 
LuRli2Si2 displays a deviation from linearity as discussed 
in Ref. 47. In this case, the initial-slope Hall coefficient 
was deduced by extrapolating the differential Hall coef- 
ficient Ra(B) = dp K {B)/dB to B = 0. Although this 
procedure yielded slightly larger values of Rn, the anal- 
ysis presented here is not affected by this offset. 

We note that the error on the absolute value of i?H 
arising from the uncertainty of the thickness of the sam- 
ples is of the order of 10%. For YbRIi2Si2, the results 
were scaled by a single factor in the temperature range 
20 K < T < 400 K to the previously published data 6 . The 
fact that this leads to a very precise match of Rn(T) in 
this temperature range accounts for the error arising from 
the thickness which enters as a factor. Consequently, the 
uncertainty of the sample thickness does not obstruct a 
detailed comparison of the different YbRli2Si2 samples. 

Measurements between 2 K and 400 K were conducted 
in a Quantum Design Physical Property Measurements 
System. For measurements down to T = 15 mK a 
3 He/ 4 He-dilution refrigerator was utilized. In this case, 
the voltages were amplified by low-temperature trans- 
formers and subsequently recorded by a standard lock-in 
technique. 



C. Experimental Results and Discussion 

1. LuRh2Si2 

The temperature dependence of the linear-response 
Hall coefficient, Rn(T), for LuRli2Si2 is depicted in 
Fig. [5] for temperatures between 17 mK and 400 K. For 
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Here, we approximated p w cr _1 which is justified given 
the small Hall angle of less than 3° . The Hall coefficient 
of the individual bands is introduced as 



J xyz 



(0 



The total resistivity pg given by 

po 1 = e p^y 1 



(ii) 



(12) 



FIG. 6: (color online) Temperature dependence of the linear- 
response Hall coefficient Rh of LuRli2Si2. The inset displays 
the resistivity of LuRri2Si2 as a function of temperature. Solid 
lines represent the simulated data according to the two-band 
model described by Eqs, 1101 to 1151 using the parameters spec- 
ified in Tab. [TT] (see text). Dashed line in the inset denotes 
the electrical resistivity calculated on the basis of a measured 
phonon spectrum^. Adding the experimentally determined 
residual resistivity pr = 1.2 \iQcm in accordance with eq. [14] 
yields a precise match with the measured data. 



LuRli2Si2 the Hall coefficient is always positive indicat- 
ing hole-like charge carriers in agreement with the pre- 
dictions of the /-core calculations on YbRli2Si2. The 
temperature dependence displays a constant value below 
20 K at i? H ~ 5.3 x 10" 10 m 3 /C followed by a crossover 
to another constant value of « 3.1 x 10 _10 m 3 /C above 
100 K. We note that a very similar behavior of Rn(T) 
has been observed for the non-magnetic La analogues of 
the CeMIns (M=Co, Ir, Rh) family of compounds^ 

For elemental copper a similar temperature depen- 
dence was observed. It was suggested to arise from two 
bands effectively contributing to the Hall coefficient, with 
their weights changing as a function of temperature^ In 
the following we demonstrate that the combination of 
band-structure calculations and comprehensive electrical 
transport measurements allow to model the temperature 
dependence of i?u quantitatively. 



was measured simultaneously with the Hall effect and is 
displayed in the inset of Fig. [5] By introducing the ratio 
r = p(l)/p(2) of the resistivities of the two bands we 
obtain the form 



i?H 



fl H (l) + r 2 fl H (2) 
(r + lf 



(13) 



Here, it becomes obvious that the overall Hall coefficient 
is only a function of the ratio r but not of the absolute val- 
ues of p(i), provided the Rn(i) are temperature indepen- 
dent. For LuRli2Si2 this latter assumption is supported 
by the band-structure calculations which yield constant 
values of Rn{i) up to 400 K, i.e., the thermal broadening 
of the Fermi surface has negligible influence, a typical 
behavior of conventional metals due to their high Fermi 
temperatures. We rather assume that merely p(i) are 
temperature dependent. For LuRli2Si2 we model the re- 
sistivity as a sum of different contributions according to 
Matthiessen's rule restricting ourself to a residual (/9r) 
and a phononic (pp) term: 



p(i) = PB,(i) + pp(i) 
The Bloch-Griineisen law 



Pp(i) 



e D /T 



sum (x) 



-dx 



(14) 



(15) 



2. Application of a two-band model to LuRh'iSi^ 

In the light of the results of our calculation (Tab. |T| 
it is reasonable to interpret the temperature dependence 
of i?H of LuRli2Si2 within a two-band model. Here, the 
observed crossover may be interpreted as the transition 
between the limits of the respective band dominating the 
total Hall coefficient. This may result from a shift of 
the relative scattering rate of the charge carriers in the 
individual bands off either phonons at high temperatures 
or static defects at low temperatures. For a quantitative 
analysis, we rewrite Eq. [3] for two bands as a function of 



describes the phononic component very well (cf. inset of 
Fig. [5]). This is corroborated by the agreement of the 
measured resistivity and the electrical resistivity calcu- 
lated using a phonon DOS derived from measured inelas- 
tic neutron scattering spectra^L (cf. Fig. [S]). In Eq. [T51 
C(i) is a constant related to the electron-phonon scatter- 
ing probability of each band, and 0d = 380 K is the De- 
bye temperature determined from specific heat^. Tak- 
ing Eqs. [TO] to [TS] together one recognizes that the to- 
tal Hall coefficient is determined at low temperatures by 
the ratio of the residual resistivities, and at high tem- 
peratures by that of the phonon scattering rates. This 
is in good agreement with the experimental data: The 
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TABLE II: Parameters calculated for LuRli2Si2 within the 
two-band model. The values were obtained from the param- 
eters of Tab. [I] following the recipe outlined in the text. In- 
serting the values listed here, we simulate Ru(T) and p(T) as 
shown by the solid lines in Fig. [6] 



i Rn(i) 




Oii) 


(10- 10 m 3 /C) 


(\iQ. cm) 


(uXi cm) 


1 21 


2.75 


7 


2 4.2 


2.11 


1.23 


resistivity ratio r 


1.3 


5.7 



low-temperature constant regime in i?n(T) is observed 
in the temperature range where the resistivity is almost 
constant. By contrast, the high-temperature regime of 
i?H (T) corresponds to a range where p(T) appears to be 
dominated by electron-phonon scattering as indicated by 
the fact that p(T) amounts to more than 10 times its 
residual value pa- Finally, the crossover is centered at 
T = 50 K where the resistivity is twice its residual value 
implying that both contributions pp and pp are equal at 
this temperature. 

Equations [TU] to [TS] contain in total six free parame- 
ters: the Hall coefficients, the residual resistivities, and 
the phonon scattering rates of the two bands. In or- 
der to fit these equations to our data we proceeded as 
follows: Firstly, we utilized the results of our band struc- 
ture calculation (Tab. Q]) in Eq. [11] to obtain the con- 
tributions Rn(i) of the individual bands. These results 
are listed in the first column of Tab. [H] Secondly, these 
Rn(i) and our experimental Rp are employed to obtain 
r from Eq. 1131 This step is performed with the value of 
i?H measured at low temperatures yielding r — 1.3, as 
well as in the high-temperature limit. However, in the 
latter case no exact solution is possible since the solu- 
tion space is limited to i?H > 3.5 x 10~ 10 m 3 /C where 
r = 5.7 for the values Rn(i) obtained on the basis of 
the calculated electronic structure. The discrepancy be- 
tween the measured and the calculated Hall coefficient 
at high temperatures may also be corrected by a change 
of R n (2) to 3.6 x 10" 10 m 3 /C This might indicate that 
the assumption of an isotropic relaxation time is not fully 
justified. However, we rather stick to the results of the 
band-structure calculation as any change would be arbi- 
trary. Thirdly, we take advantage of the fact that, at our 
lowest measurement temperature (T <C 0d), Pp is negli- 
gible leaving only pp in Eq. 1141 With r known from the 
second step and the total resistivity (Eq. [T2")) set to the ex- 
perimentally obtained value at low temperature, the indi- 
vidual pp(i) can be calculated. In the high-temperature 
regime, on the other hand, the residual term in Eq. Q3] 
is negligible and hence, the individual values of pp(i) are 
obtained from which, in turn, C(i) is inferred. All results 
are summarized in Tab. [n] 

With the parameters of Tab. [Til we are now in the po- 
sition to simulate the overall temperature dependence of 
both the Hall coefficient and the resistivity again em- 




FIG. 7: (color online) Initial-slope Hall coefficient Rh of 
YbIr2Si2- The arrow marks the temperature below which 
LFL behavior was observed^. Upper inset displays the tem- 
perature dependence of the resistivity. Lower inset sketches 
the setup. 



ploying Eqs. [TU]to[T|)] The results are included in Fig. [5] 
as solid lines. The good quantitative and the even bet- 
ter qualitative agreement of the simulated and measured 
data justify the application of the two-band model. We 
wish to emphasize that the position of the crossover in 
Rn(T) and the position where p(T) deviates from its 
residual value are not fitted but are dictated by the Debye 
temperature which was determined independently. 



3. YbIr 2 Si 2 and YbRh 2 Si 2 

For the heavy- fermion compounds YbIr2Si2 and 
YbRli2Si2 the temperature dependence of the Hall co- 
efficient is more complicated as can be seen from Figs. 
[7J and |SJ At high temperatures both compounds show a 
minimum in i?H (T) in the same temperature range where 
the resistivity assumes a maximum (cf. insets of Figs. [7] 
and[5]), namely at approximately 180 K for YbIr2Si2 and 
120 K for YbRh2Si2- This corroborates the earlier as- 
signment of this minimum in Rn(T) of YbRli2Si2 to the 
anomalous Hall effect arising from skew scattering which 
predicts such a correlation between the resistivity and 
the anomalous Hall contribution. 9 

Between 80 K and 30 K Rn(T) of YbIr2Si2 assumes a 
plateau at a value of 0.14 x 10~ 10 m 3 /C. This indicates 
that the anomalous contribution, typically being of im- 
portance around the resistivity maximum only, is super- 
posed to a normal component as expected in the theory 
of the anomalous Hall effect 4 

In the temperature range between 30 K and 8K 
a crossover to another plateau at a value of 
0.35 x 10- 10 m 3 /C is observed in R H (T) of YbIr 2 Si 2 . 
Two possible reasons may account for this observation: 
(i) The crossover might be of the same two-band nature 
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as in LuRli2Si2- However, as the crossover in YbIr2Si2 
is situated at a lower temperatures it should be accom- 
panied by a decreased value of the Debye-tenrperature. 
Unfortunately, Od is not yet known. Moreover, single 
crystals of LuIr2Si2 are not available to look for a possible 
shift of the two-band crossover in this non-magnetic ref- 
erence compound, (ii) Alternatively, the crossover might 
manifest the Fermi surface change arising from the onset 
of the Kondo screening effect which leads to itinerant / 
electrons contributing to the Fermi surface at low tem- 
peratures. 

Below 4K the Hall coefficient of YbIr 2 Si2 exhibits a 
pronounced increase, peaks at 1 K and drops at lower 
temperatures. At 0.23 K, Rn(T) changes sign and fi- 
nally saturates at the lowest temperatures at a value of 
-0.4 x 10" 10 m 3 /C. 

For YbRli2Si2 the minimum in Rn(T) at 100 K caused 
by the anomalous Hall effect is uniquely observed for all 
samples investigated. By contrast, below 50 K strong 
sample dependences are present. Figure [5] consists of 
data obtained for a large variety of samples. Three of 
theses samples were selected for low-temperature mea- 
surements down to 15 mK and represent the full range of 
sample dependences. Sample 1 and 3 exhibit a shoulder 
in Rn(T) around 15 K whereas i?n(T) of sample 2 shows 
a plateau in the temperature interval 7K < T < 20 K. 
All samples exhibit a maximum in Rn (T) around 1 K like 
in YbIr2Si2, however, at different absolute values. This 
maximum is assigned to the quantum critical spin fluctu- 
ations operating for all samples in the same temperature 
regime. Such a pronounced extremum was also reported 
for CeAflns where antiferromagnetic fluctuations were 
suggested as the microscopic origin. 48 '™ In YbRli2Si2, 
NMR investigations^ revealed antiferromagnetic fluctu- 
ations to be present in the designated temperature range. 
However, the tem per ature dependence of Rn(T) oc T _1 
predicted in Ref. [5fJ is not observed in YbRli2Si2 nor 
YbIr 2 Si 2 . 

Finally, at the lowest temperatures all samples of 
YbRli2Si2 show a saturation of i?n(T), setting in just 
below the Neel temperature. However, the saturation 
value appears to be sample dependent. 



4- Sample Dependences 

A series of YbRli2Si2 samples with different residual re- 
sistivities have been investigated above 2 K. It turns out 
that the saturation values at lowest temperature corre- 
late with the values at the plateau/ shoulder around 20 K. 
The plateau is more pronounced for samples with a lower 
saturation value: Sample 1 obeys a comparably slight 
shoulder and saturates at the highest low-temperature 
value. Sample 2 depicts the most pronounced shoulder 
resembling a plateau and exhibits the lowest saturation 
value. Sample 3 with an intermediate saturation value 
obeys a more pronounced shoulder than sample 1. The 
correlation indicates that the maximum around 1 K is 
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FIG. 8: (color online) Temperature dependence of the Hall 
coefficient for different samples. See Tab. Mil for legend. Re- 
sults for samples of the same batch are shown in identical 
color. Arrow indicates the Neel-temperature. Inset displays 
the resistivity of two selected crystals. 
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p(300K) 
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1.9 
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3.5 


X 


#3 
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3.4 


+ 




63114 


16.1 


3.7 


• 


#2 


63116 


24.0 


1.3 


> 




63116 


24.2 


1.5 


o 




63116 


26.5 


1.1 



TABLE III: Sample and batch numbers of the data sets shown 
in Fig. [8] together with the value of the Hall coefficient (in 
units of 10 _10 m 3 /C) and the resistivity ratio at 2K. 



caused by a superposed contribution which itself is not 
affected by the sample dependences. 

It is observed that samples from the same batch (cf. 
colors in Fig. 8) show almost identical Rn(T) curves. 
From this we infer that sample dependences arise from 
slight differences in crystal growth. By contrast, a corre- 
lation between sample quality and these sample depen- 
dences can not be found. This is quantitatively analyzed 
in Tab . 11 VI for the selected samples by a comparison of the 
low-temperature saturation value of i?H with the resid- 
ual resistivity ratio (RRR). An equivalent conclusion is 
found for the larger set of samples: In the absence of 
measurements in the mK-range, the resistivity ratio at 
2 K is used to quantify the sample quality which appears 
to be uncorrelated with the Hall coefficient at 2 K. 

No sample dependences were observed for LuRh2Si2 
for which three samples where investigated. In the case 
of YbIr2Si2 only one sample without indium enclosures 
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TABLE IV: Calculated and experimental Hall coefficients. 
The calculated 7?h values are derived by considering the two 
major bands with the assumption of equal relaxation times 
(see text). Values of the zero-temperature Hall coefficient, 
Rn(T — > 0), extrapolated from measurements are given for 
the related materials. _Rh is in units of (10 -10 m 3 /C). 



Calculation 




Experiment 




System Method 


Rh 


Sample 


RRR Rh 


(T-> 0) 


YbRh 2 Si 2 4/ core 


5.16 


LuRh 2 Si 2 


17 


5.3 






YbRh 2 Si 2 






YbRh 2 Si 2 4/ RBC 


-0.39 


sample 1 


70 


2.0 






sample 2 


120 


0.1 






sample 3 


40 


1.2 


YbIr 2 Si 2 4/ RBC 


-0.26 


YbIr 2 Si 2 


325 


-0.4 



could be identified. In samples with indium enclosures 
the rearrangement of the current distribution largely dis- 
turbs the Hall-effect measurement and therefore no state- 
ment on sample dependences can be made. However, the 
fact that the calculated and the measured Hall coefficient 
agree suggests that the measurements depict the intrinsic 
behavior. 



5. Comparison of Theory and Experiment 

Unfortunately, for YbIr2Si2 and YbRh2Si2 it is not 
possible to apply the two-band model as done for 
LuRli2Si2 because both the resistivity and the Hall coeffi- 
cient contain additional (quantum critical) contributions. 
Thus, the number of unknown parameters would increase 
and could in particular not be mapped with measured 
quantities. For a qualitative discussion we make the sim- 
plified assumption of equal relaxation rates for the two 
bands which yield the Hall coefficients listed in Tab. IIVI 

For YbIr2Si2, the agreement between the zero- 
temperature Hall coefficient, Rb_{T — > 0), extrapolated 
from measurements, and the calculated value is remark- 
able (see Tab. IIV|) . In the case of YbRli2Si2 our band 
structure calculations predict a value lower than the ex- 
perimental Rb_(T — > 0). This might be due to deviations 
from equal relaxation rates as the sample dependences 
indicate that small changes can have large influence. 

The most straight-forward interpretation of the sam- 
ple dependences in YbRli2Si2 arises from the insight pro- 
vided by the band-structure calculations. As shown in 
section III El the two bands dominating the Hall coeffi- 
cient are of opposite character and almost compensate 
each other. The actual value of the total Hall coefficient, 
therefore, depends sensitively on the ratio of the scatter- 
ing rates of the individual bands because they enter as a 
weighting factor in the summation of the individual con- 
tributions. Hence, it is reasonable to assign the observed 
sample dependences to changes of the relative scatter- 
ing rates. This is in agreement with the fact that other 



properties like specific heat, susceptibility and even resis- 
tivity (cf. inset of Fig. [8]) do not obey such strong sample 
dependences as none of these properties depends this sen- 
sitively on the ratio of the scattering rates. In fact, the 
resistivity is a sum of the two (Eq. [T^|) . 

Finally, the fact that samples of the same batch ex- 
hibit almost identical behavior in Rb_(T) allows us to 
surmise that the sample dependences are related to tiny 
differences in the actual stoichiometry caused by dif- 
ferent crystal growth conditions. Such sensitivity on 
minute changes of the composition is known, for instance, 
for the heavy- fermion superconductor CeCu2Si2 where it 
leads to even more dramatic effects, which include drastic 
changes in the ground stated 

The sample dependences in i?u (T) are observed to set 
in around 70 K. Below 10 K they are fully developed 
and appear to be conserved down to the lowest temper- 
atures as an offset between different samples. Conse- 
quently, the low-temperature Hall coefficient reflects the 
Fermi surface with sample-dependent, but fixed, weight 
of the individual sheets. This indicates that the Hall 
crossover, monitoring the Fermi surface reconstruction at 
the QCP, is robust against sample dependences as indeed 
observed— 

The comparison of the calculated Hall coefficient for 
the limiting cases of localized (4/ core) and itinerant (4/ 
RBC) 4/ electrons in Tab. IIVI shows that the inclusion 
of the Yb 4/ states into the Fermi volume leads to a 
decrease of the Hall coefficient. Consequently, the finding 
of a jump from larger i?u at zero field towards a lower 
value at elevated fields in isothermal scans 6,7 indicates a 
localization of the / electrons on the low field side of the 
QCP in YbRh 2 Si 2 . 



IV. CONCLUSION 

We have calculated the electronic band structure of 
YbRh 2 Si2 and YbIr 2 Si2 both with and without taking 
the Kondo scattering into account. Two bands were 
found to dominate the Hall coefficient. Both these bands 
are hole-like in the case of the /-core calculations neglect- 
ing the Kondo effect but are of opposite character for the 
case of the renormalized band calculation. The derived 
results allow for an in-depth analysis of the Hall coeffi- 
cient of the non- magnetic reference compound LuRli2Si2. 
We are able to quantitatively understand the temper- 
ature dependence of the Hall coefficient in terms of a 
two-band model. 

Furthermore, we present Hall effect measurement on 
YbIr2Si2- Here, the temperature dependence of the Hall 
coefficient parallels many features known for YbRli2Si2. 
In particular, the anomalous contribution is seen to fol- 
low the expected trend. Remarkably, the Hall coefficient 
derived from the renormalized band calculation is in very 
good agreement with the measured value at lowest tem- 
peratures. 

Finally, the sample dependences of the low- 
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temperature Hall coefficient of YbRli2Si2 are discussed in 
terms of the two bands predicted by the calculations and 
seen in LuRli2Si2. The fact that the renormalized band 
calculation predicts the two bands to almost compensate 
each other indicates that the sample dependences arise 
from small changes of the scattering rates for the individ- 
ual bands. These changes are ascribed to minute differ- 
ences in the sample composition as samples of the same 
batch show almost identical behavior. More importantly, 
despite the strong sample dependencies our comprehen- 
sive study on YbRli2Si2 confirms that the distinct change 
of the Hall coefficient, found in isothermal scans across 
the quantum critical point marks a substantial change 



of the Fermi surface in the same way as expected in the 
Kondo breakdown scenario^ 



Acknowledgments 

The authors would like to thank P. Gegenwart, S. 
Kirchncr, Q. Si, R. Valenti, T. Westerkamp and G. Wig- 
ger for fruitful discussions. This work was partially sup- 
ported by NSF-DMR-0710492, NSF-PHY-0551164, FP7- 
ERC-227378 and DFG Forschergruppe 960. 



' Electronic address: Sven.Friedemann@cpfs.mpg.de 

1 P. Gegenwart, Q. Si, and F. Steglich, Nature Phys. 4, 186 
(2008). 

2 J. Custers, P. Gegenwart, H. Wilhelm, K. Neumaier, 
Y. Tokiwa, O. Trovarelli, C. Geibel, F. Steglich, C. Pepin, 
and P. Coleman, Nature 424, 524 (2003). 

3 O. Trovarelli, C. Geibel, S. Mederle, C. Langhammer, 
F. M. Grosche, P. Gegenwart, M. Lang, G. Sparn, and 
F. Steglich, Phys. Rev. Lett. 85, 626 (2000). 

4 P. Gegenwart, J. Custers, C. Geibel, K. Neumaier, 
T. Tayama, K. Tenya, O. Trovarelli, and F. Steglich, Phys. 
Rev. Lett. 89, 056402 (2002). 

5 P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. 
Phys.: Condens. Matter 13, R723 (2001). 

6 S. Paschen, T. Liihmann, S. Wirth, P. Gegenwart, 
O. Trovarelli, C. Geibel, F. Steglich, P. Coleman, and Q. Si, 
Nature 432, 881 (2004). 

7 S. Friedemann, N. Oeschler, S. Wirth, C. Krellner, 
C. Geibel, F. Steglich, S. Paschen, S. Kirchner, and Q. Si, 
to be published (2010). 

8 A. Fert and P. M. Levy, Phys. Rev. B 36, 1907 (1987). 

9 S. Paschen, T. Liihmann, S. Wirth, O. Trovarelli, 
C. Geibel, and F. Steglich, Physica B 359-361, 44 (2005). 

10 J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 

11 A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

12 T. Moriya and T. Takimoto, J. Phys. Soc. Jpn. 64, 960 
(1995). 

13 Q. Si, M. S. Rabello, K. Ingersent, and J. L. Smith, Nature 
413, 804 (2001), cond-mat/0011477. 

14 T. Senthil, M. Vojta, and S. Sachdev, Phys. Rev. B 69, 
035111 (2004). 

15 C. Pepin, Phys. Rev. Lett. 98, 206401 (2007). 

16 G. A. Wigger, F. Baumberger, Z. X. Shen, Z. P. Yin, W. E. 
Pickett, S. Maquilon, and Z. Fisk, Phys. Rev. B 76, 035106 
(2007). 

17 M. R. Norman, Phys. Rev. B 71, 220405(R) (2005). 

18 S. Danzenbacher, Y. Kucherenko, D. V. Vyalikh, 
M. Holder, C. Laubschat, A. N. Yaresko, C. Krellner, 
Z. Hossain, C. Geibel, X. J. Zhou, et al., Phys. Rev. B 
75, 045109 (2007). 

19 T. Jeong, J. Phys.: Condens. Matter 18, 10529 (2006). 

20 S. Friedemann, N. Oeschler, C. Krellner, C. Geibel, 
S. Wirth, F. Steglich, S. Paschen, S. MaQuilon, and 
Z. Fisk, Physica B 403, 1251 (2008). 

21 U. Kohler, N. Oeschler, F. Steglich, S. Maquilon, and 



Z. Fisk, Phys. Rev. B 77, 104412 (2008). 

22 R. M. White and P. Fulde, Phys. Rev. Lett. 47, 1540 
(1981). 

23 P. Fulde and J. Jensen, Phys. Rev. B 27, 4085 (1983). 

24 G. Zwicknagl, Adv. Phys. 41, 203 (1992). 

25 P. Fulde, P. Thalmeier, and G. Zwicknagl, Solid State 
Physics, vol. 60 (Academic Press, New York, 2006). 

26 G. Zwicknagl, Phys. Scripta T49A, 34 (1993). 

27 P. Thalmeier and G. Zwicknagl, Handbook on the Physics 
and Chemistry of Rare Earths, vol. 34 (Elsevier Science 
Ltd, 2005). 

28 U. Pulst, Ph.D. thesis, TH Darmstadt (1993). 

29 O. Stockert, E. Faulhaber, G. Zwicknagl, N. Stiifier, H. S. 
Jeevan, M. Deppe, R. Borth, R. Kiichler, M. Loewenhaupt, 
C. Geibel, et al., Phys. Rev. Lett. 92, 136401 (2004). 

30 P. Thalmeier, G. Zwicknagel, O. Stockert, G. Sparn, and 
F. Steglich, in Frontiers in Superconducting Materials, 
edited by A. V. Narlikar (Springer, Berlin, Heidelberg u.a., 
2005), pp. 109-182. 

31 G. Zwicknagl, J. Low Temp. Phys. 147, 123 (2007). 

32 I. Eremin, G. Zwicknagl, P. Thalmeier, and P. Fulde, Phys. 
Rev. Lett. 101, 187001 (2008). 

33 Z. Hossain, C. Geibel, F. Weickert, T. Radu, Y. Tokiwa, 
H. Jeevan, P. Gegenwart, and F. Steglich, Phys. Rev. B 
72, 094411 (2005). 

34 O. K. Andersen, Phys. Rev. B 12, 3060 (1975). 

35 H. Skriver, The LMTO Method, vol. 41 (Springer- Verlag, 
Berlin, 1984). 

36 R. C. Albers, A. M. Boring, and N. E. Christensen, Phys. 
Rev. B 33, 8116 (1986). 

37 U. v. Barth and L. Hedin, J. Phys. C 5, 1629 (1972). 

38 C. J. Bradley and A. P. Cracknell, Mathematical Theory 
of Symmetry in Solids (Clarendon Press, Oxford, 1972). 

39 T. Jeong and W. E. Pickett, J. Phys.: Condens. Matter 
18, 6289 (2006). 

40 A. Hiess, O. Stockert, M. Koza, Z. Hossain, and C. Geibel, 
Physica B: Condensed Matter 378-380, 748 (2006). 

41 O. Stockert, M. Koza, J. Ferstl, A. Murani, C. Geibel, and 
F. Steglich, Physica B 378-380, 157 (2006). 

42 G. Zwicknagl, in Properties and Applications of Thermo- 
electric Materials, edited by V. Zlatic and A. C. Hewson 
(Springer, 2009), NATO Science for Peace and Security 
Series. 

43 P. Gegenwart, Y. Tokiwa, T. Westerkamp, F. Weickert, 
J. Custers, J. Ferstl, C. Krellner, C. Geibel, P. Kerschl, 



13 



K. H. Muller, et al., New J. Phys. 8, 171 (2006). 
P. M. C. Rourke, A. McCollam, G. Lapertot, G. Knebel, 
J. Flouquet, and S. R. Julian, Phys. Rev. Lett. 101, 237205 
(2008), 0807.3970. 

G. Knebel, R. Boursier, E. Hassinger, G. Lapertot, P. G. 
Niklowitz, A. Pourret, B. Salce, J. P. Sanchez, I. Sheikin, 
P. Bonville, et al., J. Phys. Soc. Jpn. 75, 114709 (2006). 
C. M. Hurd, The Hall effect in metals and alloys (Plenum 
Press, New York, 1972). 

S. Friedemann, N. Oeschler, S. Wirth, F. Steglich, 
S. MaQuilon, and Z. Fisk, physica status solidi 247, 723 
(2010), 0906.2649. 

M. F. Hundley, A. Malinowski, P. G. Pagliuso, J. L. Sarrao, 
and J. D. Thompson, Phys. Rev. B 70, 035113 (2004). 



J. Ferstl, Ph.D. thesis, Technische Universitat Dresden 
(2007). 

Y. Nakajima, K. Izawa, Y. Matsuda, K. Behnia, H. Kon- 
tani, M. Hedo, Y. Uwatoko, T. Matsumoto, H. Shishido, 
R. Settai, et al., J. Phys. Soc. Jpn. 75, 023705 (2006). 
K. Ishida, K. Okamoto, Y. Kawasaki, Y. Kitaoka, 
O. Trovarelli, C. Geibel, and F. Steglich, Phys. Rev. Lett. 
89, 107202 (2002). 

F. Steglich, P. Gegenwart, C. Geibel, P. Hinze, M. Lang, 
C. Langhammer, G. Sparn, T. Tayama, O. Trovarelli, 
N. Sato, et al., in More is different: Fifty years of con- 
densed matter physics, edited by N. P. Ong and R. N. 
Bhatt (Princeton University Press, Princeton, 2001), pp. 
191-210. 



